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Most previous theoretical studies on the optical torque exerted by light on dipolar particles are 
incomplete. Here we establish the equations for the time-averaged optical torque on dipolar bi¬ 
isotropic particles. Due to the interference of scattered fields, it has a term additional to that 
commonly employed in theory and experiments. Its consequences for conservation of energy, angular 
momentum, and effects like negative torques, are discussed. 


In recent times, and with the steady improvement of 
particle manipulation techniques, the torque exerted by 
light on particles is gaining increasing attention. How¬ 
ever, in spite of the last advances attained on the theory 
and experiment interpretation of the optical force (OF) 
on small particles [iMtI, there is no yet aparallel one on 
the optical torque (OT). Most studies [8|,l9j have no quan¬ 
titative assesment with theory or, if theoretical, involve 
a static formulation [^, ITfil - fl^ and thus are incomplete 
and hence incorrect: they do not comply neither with 
the conservation of angular momentum nor with that of 
energy. 


We know few exceptions: reports involving Mie’s calcu¬ 
lations for larger particles or sets of spheres [lllll . How¬ 
ever, the latter convey computational procedures that, al¬ 
though yielding rigorous exact resuts, do not completely 
allow to see through analytic calculations the details of 
the different involved physical processes. These, like in 
the case of optical forces, are clearly illustrated in the 
equations for the torque on a dipolar object that, being 
a canonical system, is of great utility both for forthcom¬ 
ing incresingly controlled experiments and interpretation 
purposes. This has already been extensively seen in the 
study of photonic forces. 


In this letter we establish the theory of the OT on dipo¬ 
lar particles. By ’’dipolar” we mean those in the wide 
sense, i.e. whose polarizabiliities are fully given by the 
first electric and magnetic partial waves [7|,ll5|, (e.g. Mie 
coefficients if they are spheres). We put forward that the 
time-averaged OT < F > on a dipolar particle contains 
a component additional to that < Fq > (named intrin¬ 
sic torque) commonly derived from static equations and 
that, for reasons later discussed, coincides with what we 
shall name extinction optical torque (EOT). This addi¬ 
tional component is the recoil optical torque (ROT) due 
to interference of fields scattered by the particle, and is 
essential to fulfill the conservation of angular momentum 
and energy. It is curious, on the other hand, that this 
ROT implicitely appears as a well-known textbook re¬ 
sult 0 for the rate of angular momentum radiated by 
a dipole, even though apparently without being noticed 
by many modern research works. However, the ROT is 
essential (as noticed in 0 for sets of large particles) 


to predict new OT effects, such as e.g. the appearence 
of negative optical torques (NOT). A phenomenon which 
keeps an analogy with pulling forces [ITHlSj, and that it 
is much more common than the latter, nor requiring a 
special combination of incident beam structure and ob¬ 
ject composition and morphology. 

We consider monochromatic light fields, so that 
the electric and magnetic vectors are: £ (r, t) and 
B{r,t): £{r,t) = 5J[E(r) exp(—iwt)] and B(v,t) = 

5J[B(r) exp(—iwt)]. K denotes real part. We will ad¬ 
dress time-averages of electromagnetic quantities, in¬ 
cluding the electric and magnetic energy densities, re¬ 
spectively: < We >= (e/I67r)|Ep and < Wm >= 
(l/167ru)|BP, and the time-averaged Poynting vector: 

< S >= (c/87r^)K(E X B*). 

Let a particle of index Up be embedded in a medium 
with dielectric permittivity e and magnetic permeability 
pL, illuminated by an electromagnetic field whose elec¬ 
tric and magnetic vectors are: E^*), The opti¬ 

cal torque F = dJ/dt, exerted by this field on the ob¬ 
ject, comes from the conservation of the angular momen¬ 
tum J = r X S , r denoting a position vector of space. 
On time averaging, such conservation law: < F >=< 

dJrnech/dt > -|- < d3fieid/dt >= f^rx < T > els' 0 

leads to the torque time-average: 

< F >= —/ r X [e(E • s)E* 

Stt Js 

+M-\B • s)B* - i(e|E|2 + ^-i|B|2)s]d5. (r = rs). (1) 

The angular momentum < J > has been splitted into 
the mechanical one: < Jmech > and that of the wave: 

< J field >, (< d3 field/dt >= 0 for these fields). T is 
Maxwell’s stress tensor. dS denotes the element of any 
surface S (which will be considered a sphere for calcu¬ 
lation purposes) that encloses the particle, and s is its 
local unit outward normal, r = rs is a point of S. 

The wavefields in Eq.(IT|) are total fields: the sum of 
those incident and scattered, denoted by the {i) and (s)- 
superscripts: E^*) -l-E^'’^ B^'^ -l-B^'’^. Obviously the last 
two terms of the integrand of o do not contribute to 

< F >. 

Let the particle be magnetodielectric and bi-isotropic 
[13, dipolar in the wide sense, for example if a sphere [3, 
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Il5l| , its magnetodielectric response is characterized by its 
electric, magnetic, and magnetoelectric polarizabilities : 

^ 2k^ ^ 2k^ ’ ^em ^ ’ ^me ^ ' 

ai, 5i and Ci standing for the electric, magnetic, and 
magnetoelectric 2l| first Mie coefficients, respectively. 
Then the electric and magnetic dipole moments, p and 
m, induced on the particle by the illuminating field, are 




Bd) 


= ( 2 ) 


First consider an incident plane wave, circularly po¬ 
larized (CP), (cf. Fig.l): Ed) = e,e*'=d. r)^ Bd) = 
b^g*fc(si-r)^ The scattered fields at a point P: rs are 


1 1 7 A’ 

= -[3s(s-p)-p](^- ^)exp(iA:r) 


6 r 


( 3 ) 


1 ih 

= /r[3s(s • m) — m](—-^) exp(zA:r) 

&„xp)£:;:+b(,,)5^E<ihi, . 

V e r 

With the far-zone scattering amplitudes being 


( 4 ) 


e{sxy) = *^[-(s X p) X s - \j-{s X m)], (5) 


b(sa;y) = X m) X s-I-X p)]. (6) 

s = (sxy, Sz), bi = nSi XBi, Bi-Si = hi-Si = 0; b = ns x e, 
e • s = b • s = 0. The rest of terms of ([3]) and (jl]) are the 
fields in the near and intermediate zones 161. 


X 



FIG. 1: (Color online). Scattering geometry illustrating a 
negative optical torque of a CP plane wave on two particles. 

We adopt a Cartesian framework, (cf. Fig.l and con¬ 
sider by now only one dipolar particle in 0), where the 


incident wave has = (0,0,1), and = e(l,±i,0), 
bi = ne(=Fi,l,0); n = y/eJI. The upper and lower sign 
stands for left (LCP) and right (RCP) circular polar¬ 
ization. Introducing Eqs. © - © into o, writing 
dS = sin OdOdcj), the interference between the incident 
and scattered fields in the first two terms of m yield in¬ 
tegrands: [(p* X s)(ei • s) -I- 2(p- s)(s X e*) -t- (m* x s)(bi • 
s) -I- 2(m • s)(s X b*)] which come from the near-zone 
scattered fields in the products: r^[e(s x e)®)*)(E)®) -s) -|- 
e(s X E«*)(E)®) • s) -k n-\s X B)®)*)(BW • s) -k ^-^(s x 
B(*)*)(bC) . s)]. Notice that these r~^ near-zone scat¬ 
tered fields [cf. Eqs. © and ©] are those for which 
there is cancellation of the factor in the integrand of 
o, and hence lead to results independent of the radius 
of integration r. By contrast, the contribution of far and 
middle-zone terms give results proportional to and r, , 
which are zero allowing the integration sphere to be arbi¬ 
trarily close to the enclosed particle, which is treated like 
a point. Thus the integrals over the polar and azimuthal 
angles lead to 

< To >= ijJKp X eW*) -k (m X B«*)]. (7) 

< Eq > arising from the interference of the incident and 
near-zone scattered fields (in particular, when extremely 
close to the object they are static) comes from the ex¬ 
tinction of the incident angular momentum, (in fact of 
the spin part, SAM), and its transference to the particle. 
For an incident plane wave there is no contribution of 
the orbital angular momentum (0AM) to this extinction 
optical torque (EOT) < Eg >. 

On the other hand, the interference of the scattered 
fields in the integrand of © contain terms: 2[(p-s)(p*x 
s) -I- (m • s)(m* X s)] which arise from the product 
of the middle-zone with the far-zone fields with an r- 
dependence: r~^ andr“^, respectively, [cf. Eqs. © - © 
or © - ©]: r3[e(sxE(®)*)(E)®)-s)-kM-HsxB(®)*)(B(®)- 
s), again cancelling the factor of the integrand. After 
angular integration, they lead to the scattering, or recoil, 
contribution to the optical torque: 

< E^®^ >= —■^[^^(p* X p) -I- ^3(m* X m)]. (8) 

O 6 

5 means imaginary part. As said, < E^®) > equals the 
rate of radiation of electromagnetic angular momentum 
to infinity < dJ/dt > from the oscillating electric and 
magnetic dipoles 0- This is logical on taking into ac¬ 
count that these dipoles are induced on the particle by 
the incident wave and then re-radiate. To this compo¬ 
nent, half is contributed by the 0AM and the other half 
by the SAM. 

Therefore, in summary, o, © and © yield the total 
time-averaged torque on the particle due to an incident 
plane wave as 

< E > = < Eg > -f < E(®) > . 


(9) 
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Which contains the contribution of the scattered, or re¬ 
coil, component < >, Eq.®, added to < Fq >, 

Eq.(I71). The former, like the latter, is an intrinsic torque. 

Eq. m is the first main contribution of this letter. 

In this regard it is interesting that the ROT ([5]) keeps 
a formal analogy with the electric-magnetic component 
[given by 5i(p x m*)] of the electromagnetic force on a 
magnetodielectric dipolar particl^as it also arises by in- 


Ctk 

terference of the scattered fields [Tj. On the other hand, 
the term < Fq >, Eq. 0 , is formally analogous to the 
time-averaged Lorentz’s force on an electric (magnetic) 
dipole from a time-varying magnetic (electric) field, given 
in terms of S[p x b*], (3[m x e*]). Thus the RT plays a 
role in the OT as important as the e—m electric-magnetic 
interaction force on an electric-magnetic dipole system. 
As a matter of fact, like the omission of the e — m inter¬ 
action avoids the prediction of phenomena like Kerker’s 
effects or pulling forces 17 that of the RT prevents 
the obtention of a NOT, i.e. a torque on the body coun¬ 
terclockwise with respect to the helicity of the incident 
light. 

Let us assume the particle being chiral, then Ogm = 
—ame- Introducing now Eqs. m into 0 - ® and using 
the optical theorem, which expressing the conservation 
of energy imposes a well-known condition (cf. Eq.(16) of 
i) between polarizabilities, which for a chiral particle 
generalizes to: 


e -—- H + 2/r|am,eP T 

47rfe 3 

2n^[e~^alame + fJ-amea*^] 

= -h {ai + (10) 

representing the particle absorption cross-section), 
we write the torque on the chiral paticle due to a CP 
incident plane wave, as: 


< F >= ±ze e 


cr(“) 

47rfc 


( 11 ) 


where z is the unit vector along OZ and the signs plus 
and minus correspond to LCP and RCP, respectively. 

Hence the OT on one bi-isotropic (or in particular, chi¬ 
ral) body illuminated by a plane wave is essentially mea¬ 
sured by its absorption cross-section. This coincides with 
previous knowled ge d erived both from Mie theory and 
experiments [^, lol Il3|| in large purely dielectric isotropic 
spheres. Therefore, (ED shows that even if the parti¬ 
cle is magnetodielectric and bi-isotropic, being dipolar, 
it experiences an OT due to a plane wave which is es¬ 
sentially given by the absorption cross section, thus the 
OT being zero if the particle is non absorbing. Unlike 
the OF, the OT does not change depending on whether 
the particle satisfies Kerker’s |22|,[23j or any other magne¬ 
todielectric conditions. Only the absorption cross section 
plays a role. As seen, this characteristic is a consequence 
of the conservation of both energy, ruling Eq. (nnD, and 


angular momentum, leading to O- This is the second 
contribution of this letter. 

The aforementioned duality of the OT torque and OF 
equations on dipolar particles is even more notorious 
when the incident field is an arbitrary wave. To do this, it 
is convenient to represent the fields, both in the near and 
intermediate regions, by the angular spectrum of plane 
waves, 


24? 


In this case a calculation with 0 shows that Eq. (|9D 
applies with the ROT < F*^®^ > remaining as in (jSD, how¬ 
ever the EOT now becomes instead of < Fq >: 


< f(°) >=< Fo > 

+ ^y^3{^(p • V)BW* - /r(m • V)eW*}. (12) 

Where now the contribution from the 0AM to the second 
term of m is half that from the SAM, and < Fq >, con¬ 
tributed by the SAM, is given by Eq. 0 with the vectors 
and of course being those of the arbitrary inci¬ 
dent field. Notice that the additional terms appearing in 
m are the formally analogous to those of the ’’dipolar 
force” which in the OF add to the ’’Lorentz’s force terms, 
and are due to the spatial structure of the incident field. 
This is further seen by writing (ED as 

< r'»> >= 

2 

(13) 

Eqs. (TT^ and (1131) are the third main contribution of 
this letter. 

Eq.(ED has a formal analogy with the well-known ex¬ 
pression for the extinction OF (i.e. that containing the 
electric and the magnetic OF) on a magnetodielectric 
dipolar object 0. Indeed, adding (fT^ and 0 in 0, 
the analogy of the OT and the OF on a magnetodielec¬ 
tric dipolar particle is complete as long as the duality of 
factors p and m, and E^*) and B(®\ appearing in the OT 
and in the OF is concerned. It should be stressed that 
now < F > is no longer proportional to 

Next, we illustrate this analysis addressing the phe¬ 
nomenon of NOT in sets of dipolar particles. We con¬ 
sider two identical dielectric spheres with centers at 
ro = (0,0,0) and rg = (0, R, 0), (cf. Fig. I). Illumi¬ 
nation occurs with a plane CP wave. Considering the 
coupling between the dipoles induced in both particles, 
the resulting OT on the two sphere system is: 

< F >= ±ze"|ae|"R{(^P|ae|2 + e^)R + afl}. (14) 
3 47rfc 

The upper and lower sign of ± applying to LCP and RCP, 
respectively. R and I denote real and imaginary parts. 

Or, by using (ED 

< F >= ±ze^\ae\^n{ain -f afl}. 


(15) 
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Where 

cy^ikR c2^AkR 

vw = ^-L = 

are the tensor Green function eigenvalues, so that the real 
and imaginary parts of 77|| — 77 _l are TZ = 5i[?7|| — = 

|^fcii/||(fci?) + ■^{2kRsinkR + ZcoskR). X = 3 [ 77 || — 
ryl] = ^[/||(fci?) - fdkR)] - 



FIG. 2: (Color online). Averaged optical torque, normalized 
to incident energy flow < |Si| >, versus kR on two Si particles 
{Up = 3.5) separated a distance R. {a^ = 7.91 x 10® nm^, 
ai = 1.1 X lO^nm®, c = 3 x lO^^nm • s“^). The particles 
rotate with respect to an axis parallel to the OZ direction, at 
y = R/2. There are regions where < F > is negative. 

It is curious that /|| (x) and /j_ (x) are the parts (par¬ 
allel and perpendicular, respectively, to the the OY seg¬ 
ment joining both particles) of the diagonal elements of 
the correlation tensor fij(x) (bJ = 1)2,3 = x,y,z) of 
an statistically homogeneous and isotropic random field 
E’'(r), namely fuix) =< > / < \El{r)\'^ >, 

X = |r - r'|. So^at /||(x) = fyy{x) and /_l(x) = 
fxxix) = fzzix) [ii: /||(x) = 3(smx - Xcosx)/x^; 
/j_(x) = (3/2) [sinx/x - {sinx - Xcosx)/x\ Obvi¬ 
ously, when the particles approach each other, one has 
7^ ^ 2,/eR^, 1^0, {kR 0). 

One may wonder on the appearence of fijix) in this 
theory, however, this is not surprising as these functions 
are associated to t^h and rj±, and appear in any interaction 
with a kR dependence characterized by coupled dipoles. 
Notice, on the other hand, the oscillating behavior of 
the OT exerted on the two spheres versus kR, according 
to Eqs. da or da; this encompasses the existence of a 
NOT on the system. 

As an example, we consider a LCP plane wave incident 
on two Si spheres of radius 230nm in air at A ~ lS50nm, 
where the electric dipole dominates [l^. As seen there 
are separations R leading to a negative OT. of influences 
mucho more than the results in (1151) . This torque (cf. 
Fig. 2) is enhanced in the electric dipole resonance region 


in this kind of high Up particles since, on comparison, we 
obtained a torque one order of magnitude smaller in the 
equivalent region (A = TOOnm) if being of the same size, 
they are both of latex, (rip = 1.5). 

In conclusion, the expression of the ROT is general for 
any illuminating field. Our generalized EOT, together 
with the ROT, may give rise for specially structured in¬ 
cident beams, to negative torques on one dipolar object. 
Further studies for other sets of magnetodielectric, bi¬ 
isotropic, coupled dipoles remain to be done. 
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